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Abstract 

A new approach for derivation of Benney-like momentum chains and integrable 
hydrodynamic type systems is presented. New integrable hydrodynamic chains are 
constructed, all their reductions are described and integrated. New (2+1) integrable 
hydrodynamic type systems are found. 
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1 Introduction 

The integrable hydrodynamic chain 

dtAk = d^Ak+i + kAk-iAo^^, /c = 0, 1, ... (1) 

for the first time was introduced by D.J. Benney in a theory of finite-depth fluid (see p]). 
Here moments Ak are infinite number of field variables. Later, it was shown that these 
moments satisfy a dispersionless limit of KP hierarchy determined by the Sato pseudo- 
differential operator 

which in dispersionless limit is reduced to 

A = /i+^ + ^ + ... (2) 

The Benney momentum chain can be written in equivalent form (see |18j ) 

dX u?' 

h - f^K = ^[/^t - Y + ^o)]. (3) 

If A = const, then is a generating function (with respect to the parameter A) of the 
conservation law densities 

f^, = d,{^ + Ao). (4) 



In the case, when several first moments are functionally-independent {k = 0, 1, A^ — 
1), the corresponding hydrodynamic reductions ("hydrodynamic" reduction means that 
all higher moments [k = 1, 2, ...) could not depend of any derivatives of lower moments 
Ak {k = 0, 1, — 1); in opposite case, such reductions one can call as "differential" 

reductions. In this article we concentrate attention at hydrodynamic reductions only in 
spirit of pni) are the hydrodynamic type systems written in Riemann invariants r* 

rj =/i,(rK, 2 = 1, 2, AT, (5) 

i.e. hydrodynamic type system has diagonal form in these field variables and there is no 
summation over each repeated index, see for instance |HE]- The Riemann invariants 
and the characteristic velocities /ij(r) are determined by conditions 

r=H:-\ \ 5- H T + 1 = ^ + 2^ + 3^ + ... 

f^i K K K K K 

(see (j2I) and Q). These hydrodynamic type systems are integrable too (all moments A^. 
are some functions of r\ which are determined by compatibility conditions with whole 
Benney momentum chain, see |19 ). These hydrodynamic type systems have the same 
generating functions of conservation laws (see (jH)) and the commuting flows (see ^Tj; 
"commuting flows" means, that the Riemann invariants r* simultaneously are functions 
of infinite number of independent variables t^, k = 0, 1, ... , here to = x, ti = t.) 

/^(^)r(A) = 9^ ln[/x(A) - /i(A)] , (6) 

where ^ ^ 

^r(A) = 9to + jdt, + ^dt, + ... 
A A 

Moreover, the generating function (with respect to the parameter A) of solutions for any 
reduction ^ can be found by the Tsarev generalized hodograph method ( 36j, also see 

m) ^ 

x + /i,(r)t = — t = l,2,...,N (7) 

Thus, if some hydrodynamic type system is recognized as a reduction of the Benney 
momentum chain, it means that this system has most properties of the Benney momentum 
chain. 

The idea presented in this paper is the following: if one can introduce the moments 
for given integrable hydrodynamic type system (jSJ, then one can ignore the origin 
(i.e. given hydrodynamic type system) of this hydrodynamic chain 

dtA = F{A)A^, 

where A is an infinite- number component vector, F{A) is an infinite- number component 
matrix. The next step is a description of all possible integrable hydrodynamic reductions 
(one of them, of course, must be the original hydrodynamic type system 

rl = Viiv)rl, 1 = 1,2, M, (8) 
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where M is not connected with N (see ©), and Vi satisfy some nonhnear system of 
PDE's (see below). Also, we assume that Vi ^ Vk for any i ^ k (this is a necessary 
condition for the application of the Tsarev generalized hodograph method). Thus, every 
hydrodynamic chain constructed in this way can be regarded as a huge box for some 
variety of the integrable hydrodynamic type systems. 
For instance, the particular case of gas dynamics 

ut = d^[— + -], r]t = d^{ur]), (9) 

/ 7 — i 

for 7 = 2 (shallow water equations): 

= + r?], Vt = dx{uri) (10) 

satisfies for Benney momentum chain (PJ) if one introduces the moments Ak = u^rj. 

It is easy to check, that the Benney momentum chain has a more general (the Za- 

N 

kharov) reduction Ak = J2'^iVi (see [HH]), which create a dispersionless limit of vector 

i=l 

nonlinear Shrodinger equation (VNLS) 

9 N 

dtUi = dx[-^ + 22Vk], dtr]i = dx{uir]i), i = l,2,...,N (11) 

k=l 



Remark 1 Infinite series ^ under this Zakharov reduction yields a more compact ex- 
pression (see lEEl) 

N 

A = /^ + E^- (12) 

It is easy to check that the dispersionless limit of VNLS satisfies equation (0) with respect 
to equation of Riemann surface / lijj)) . 

Obviously, in both above-mentioned cases, corresponding hydrodynamic type systems 
()10|) and (fTTj) have the same generating functions of conservation law densities (j3)) and 
commuting flows (jH)) as whole Benney momentum chain (^Q). Here we demonstrate this 
approach on an example of a new hydrodynamic chain, which contains some important 
reductions well known in mathematics, fluid dynamics, nonlinear optics, biology and 
chemistry. 

The main classification problem in the theory of integrable hydrodynamic type systems 
can be re-formulated as the problem of description of all possible integrable hydrodynamic 
chains. For the simplicity, any hierarchy of the hydrodynamic chains can be written in a 
conservative form (see, for instance, ()30p below) 

dt^Ak = d^FkiAk+n, Ak+n-i, ... , Aq), k,n = 0, 1, ... 
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If one can classify all possible functions F^, it would mean that all the hydrodynamic type 
systems embedded in such hydrodynamic chains by different reductions are classified too. 
In simplest case (A^ = 2) 

dtAo = d,Fo{Ao, Ai), dtAi = d,Fi{Ao, A^, A2), 

where A2{Aq, Ai) is conservation law density of reduced hydrodynamic type system. Then 
governing equation for function A2 is quasilinear 

fvWuu = [fu - V^w - <^w^v]Wuv + [<^u + y^wWulWyv, 

where 

u = Aq, V = Ai, w = A2, f = Fo, ip = Fi. 

When Fq = v, Fi = w{u, v) — v? /2, this is 2— component reduction of the Benney 
momentum chain (see [T?l] ) 

Wuu = -WvWuv + {Wu - u)Wvv] 

if Fo = — V?, Fi = w{u, v) — uv, then corresponding equation 

Wuu = -{u + Wy)Wuv + {Wu - v)Wvv, 

can be solved in parametric form 

w = l[A"{s) + B"{r)f + [A"{s) + B"{r)][A'{s) - sA"{s) + B'{r) - rB"{r)] + 
- 2sA'{s) + 2A{s) + r'^B"{r) - 2rB'{r) + 25(r), 

V = \[A'\s) + B"{r)f + A\s) - sA"{s) + B'{r) - rB"{r), u = A"{s) + B"(r), 

where A{s) and B{r) are arbitrary functions. Thus, 2— component reduced (hydrodynamic 
tj^e) system in Riemann invariants is 

n = (A"(s) + B"{r) + r)r,, St = {A'{s) + B"{r) + s)s^. 

This system is natural 2— parametric generalization of gas dynamics (see below). 

The goal of this paper is a complete description of A^— component generalization of 
above-mentioned formulas. 

In section 2 of this paper, so-called "e— systems" are introduced. All their properties 
like conservation laws and commuting flows are described. The corresponding hydrody- 
namic chain is found by the natural introduction of moments. 

In section 3, some properties such transformations between different representations 
of this hydrodynamic chain are obtained. 

In section 4, all possible hydrodynamic reductions are found. Particular and important 
reductions of this hydrodynamic chain are emphasized. 

In section 5, a generating functions of conservation law densities, commuting flows 
and solutions (by the Tsarev generalized hodograph method) for these hydrodynamic 
reductions are constructed. 

In section 6, a general solution of these hydrodynamic type systems is presented. 
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In section 7, new (2 + 1) integrable hydro dynamic type systems are found. 
In section 8, another hydrodynamic chain is presented, all its hydrodynamic reductions 
are described. 

In section 9, we discuss some still open problems: Hamiltonian structures and inte- 
grable dispersive extensions of hydrodynamic chains and their reductions. 

In section 10 (Conclusion), we describe a general situation in theory of hydrodynamic 
chains. 

2 "e— systems". 

This class of integrable hydrodynamic type systems 

N 

ri = [r'-eJ2r^]ri, z = l,2,...,N, (13) 

m=l 

where e is an arbitrary constant, was established in |^ (also see [TH], [2S1, IHII], [SSI)- 
These hydrodynamic type systems (fT^ and its commuting flows (see below) we shall call 
"e— systems". The particular case N = 2 plays important role in gas dynamics (see Q, 
where the adiabatic index 7 = frff)) in field theory (7 = 1, Born-Infeld equation), in 
nonlinear optics (7 = 2, the dispersionless limit of nonlinear Shrodinger equation, see Q) 
and in fluid dynamics (7 = 4, the dispersionless limit of the second commuting flow to the 
Boussinesq equation). Also, "e— systems" (for arbitrary A^) are well known in differential 
geometry {e = —1/2, elliptic coordinates, see, for instance, dispersionless limit of 
Coupled KdV, see, for instance, [13]), in soliton theory {e = 1, some particular solutions 
of linearly-degenerated systems are multi-gap solutions of KdV, see [IS]), in biology and 
chemistry (e = — 1, chromatography, electrophoresis, isotahophoresis). Moreover, a gen- 
eral solution can be found explicitly (see ^2^), for instance, in one-atomic (7 = |, e = —1), 
two-atomic (7 = |, e = —2) gases (see ()13p ) and their generalization for arbitrary and 
arbitrary integer e. Thus, obvious aim is to extend a class of integrable hydrodynamic 
type systems starting from (fT^ with preservation of some properties. 

The hydrodynamic type system (|13|) has a generation function /z of conservation law 
densities; when \ 00 

N 

^ = JJ (1 _ Tm/X)-' = 1 + ai/X + as/A' + (14) 

m=l 

when A ^ (up to constant multiplier) 

N 

/i= J](r„-A)-^ = 6o + feiA + 62A' + ... (15) 

m=l 

The first series ()14p is a series of polynomial conservation law densities with respect 
to Riemann invariants (this is analogue of Kruskal series of conservation law densities 
for integrable dispersive systems like Korteweg de Vries equation). We shall call them 
as "higher" (or "positive") conservation law densities correspondingly their homogeneity. 
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Coefficients bk we shall call as "lower" (or "negative") conservation law densities (they play 
role as a " new" conservation law densities appearing under Miura type transformation in 
theory of integrable dispersive systems). 

It is easy to check, that any commuting flow (so, every Riemann invariant r* is a 
function of three independent variables x, t, r) 

rl = wl)ir)rl (16) 

to hydrodynamic type system (fT^ has velocities 

^,)(r)=9./.(_,), (17) 

where /i(-e) is some conservation law density of "(—£:)— system" 



N 



7 r 

r, = r 



+ eJ2nrl ^ = 1, 2, iV. 



m=l 



Since, "e— systems" and " (— £)-systems" have generating functions of conservation law 
densities such that 

/^(^) ■ = 1, (18) 

then a generating function of commuting flows to ()13p in Riemann invariants is (see ()16|) 
and (IT7j) l 

r\i, = ^—^ri, (19) 

(l-rVA)/i " 

and in conservative form (sf. ^ is 

/^.(A) = T^Mj)^ (20) 

where fl = /i(A) (see (HH) and (US))). 

Higher commuting flows can be obtained from (see ()14|) ) 

AT 

/i = n ~ r„/A)-" = 1 + ai/A + aa/A^ + ... 

m=l 

and formal series ^ ^ 

^r(A) = dto + ^dt, + ^dt, + ... 

A 

when A ^ oo. The corresponding generating functions of conservation laws are 

k 

a4,/i = a,[/xJ]fi„A'=-"^], A; = 0,1,2,... (21) 

m=0 

where 

n-l 

5o = C^O = I5 ^1 = 0.n = ~(^n ~ ^^^Qm-Qn-A:; n = 2, 3, ... 

?Tl=l 
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The corresponding higher commuting flows (in Riemann invariants) are 

k 

= [Y^aUr'f-nrl k = 0,l, 2, ... (22) 

m=0 

The lower commuting flows can be obtained from (see ()15|)) 

N 

/i = JJ (r„ - A)-" = bo + h~X + + ... 



~ 2 ~ 3 

d^r^) = Xdt_, + X dt_, + A (9t„3 + 



m=l 

and formal series 

r(A) 

when A — ^ 0. The corresponding generating functions of conservation laws are 

k 

dt.,_,fi = d.ifiY.bmX'^''"^ k = 0,l, 2, ... , (23) 

m=0 

where 

1 . 1 ''-K 

bo = 7-, bk = —— y bmbk-m, k = 1, 2, ... 
bo 

m=0 

The corresponding lower commuting flows (in Riemann invariants) are 

k 

dt.,_y = [Y}m{rT-'-'Vl. = 0, 1, 2, ... , (24) 

m=0 

If A ^ 00, A — >■ cxD, all the higher conservation laws for the higher commuting flows are 

k 

dt^am = d^y^asak+m-s], k = 0, 1, 2, ... (25) 

s=0 

If A — > 0, A — cxD, all the lower conservation laws for the higher commuting flows are 

k n 

dt„bk = d^^^b,an+s-k], k ^ n, dtjk = d^[^dsbk+s-n], k^n. (26) 

s=0 s=0 

If A — > 0, A — s> 0, all the lower conservation laws for lower commuting flows are 

n 

dt-r.-ibm = d^\^bkbn+m+l-k], u = 0, 1, 2, ... (27) 
A:=0 

If A ^ oo, A ^ 0, all the higher conservation laws for the lower commuting flows are 

m n 

dt_„_iam+i = d^[^aj)n+s-m], m^n, 5t_„_,a„+i = d^[^bsam+s-n], m ^ n. (28) 

fe=0 fc=0 



All these above formulas can be easily checked by a direct calculations. 

For instance, the initial system (fTSjl has the generating function of conservation laws 

dtfi = d4{\-ai)fx\, (29) 

where an infinite set of the positive (polynomial) conservation laws is 

dtittk = <9^K+i - «i«fc], k = l,2, ... (30) 

and an infinite set of the negative conservation laws is 

dt^bo = d^{-aibo), dtj)k = d^[bk^i - aibk], k = 1, 2, ... (31) 

The second commuting flow (see (fT^ and ((221)) 



■i2 



N 2 ^ ^ 



= [{r^ - er^Y.'- + -(J^r™)^ - -Y^inY.. ^ = 1, 2, ... , N 



' 2- 

m=0 1 

has the generating function of conservation laws (see (HU) 

^^t2 = dx[{>^'^ - fliA + al- a2)/i], 
where an infinite set of the positive (polynomial) conservation laws is (see (j25|) ) 

dt^ttk = d^[ak+2 - aictk+i + (a? - ct2)afc]- (32) 
an infinite set of the negative conservation laws is (see (j26p ) 

dt^bo = c?x[&o(ai - ^2)], dp,bi = d^biiaj - 02) - aibo], 
dtibk = d^[bk{al - - aibk^i + bk^2\, k = 2,3,... 

The first negative (A 0) commuting flow is (see (|24j) ) 

N 

1 J]/''"^' 

where the generating function of conservation laws is (see 



= dx^, (34) 



where an infinite set of negative conservation laws is (see ()27|) ) 

dt.,bk = 9,^, = 0, 1, 2, ... (35) 

bo 

and an infinite set of positive conservation laws is (see (j2HI)) 

dt_^ai = dr,^, dt_^ak+i = dx'^, A; = 1, 2, ... (36) 

Oq Oq 



Remark 2 Reciprocal transformation (see the first equation in 0^1}) ) 

dy_i = h^dx — aih^dti, dz = dti, 

connects system / f73')) (also, see / f37^) and / f37)) ) and its commuting flow (also, see / f3'3)) 
and ^)). 

Reciprocal transformation (see ^EW) 

dy = fi[dx + (A — ai)dti], dz = dti 

connects system i f73|) and generating function of its commuting flows (see /f7 ^) -/f ^) ) 

dtai = -A9^(l//i). 

3 New hydrodynamic chain 

The hydrodynamic type system (fT^ can be rewritten as an infinite momentum chain 

dtCk = d^Ck+i - cid^Ck, /c = 0, ± 1, ± 2, ... (37) 
where N first moments Cfc {k = 1, 2, A^) are functionally independent 

N N 

co = eJ2^nr^, = |5^(r'-)^ A; = ±l, ±2,... (38) 

m=l m=l 

Thus, all a/s, 6^ and Ck can be expressed via each other (see below). This is invertible 
transformation. 

However, now we can start our investigation namely from hydrodynamic chain written 
in form (j37|) or, for instance, (j3(Jj) without reference on original hydrodynamic type system 
()13|) (and explicit expressions (jHHjl ). If now we restrict our infinite momentum chain to 
iV— component case, then just one particular solution obviously is the hydrodynamic type 
system (fT^ (also, see (IHHjl). How to find all other possible reductions? The answer will 
be done in next section. 

Remark 3 For the first time, this hydrodynamic chain i3(J\} has been derived (in an- 
other terms) by S.J. Alber (see f^; also it has been independently obtained in another 
context by V.G. Mikhalev, see 1261), and recently by L.M. Alonso and A.B. Shabat (see 
they describe mostly differential reductions and very particular hydrodynamic re- 
ductions, in our article we describe all possible hydrodynamic reductions). Actually, 
the starting point of their investigations was the hydrodynamic type "e— system" , when 
e = —1/2, related by generalized reciprocal transformation with averaged (by the Whitham 
method) integrable systems (determined by scalar second order spectral transform with 
energy- dependent potential, see also mentioned references) related with hyperelliptic sur- 
faces. 
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The generating function of conservation laws for the hydrodynamic chain (|37|) is ex- 
actly where (sf. jH))) 

oo oo 

/i = 1 + ^OfcA"'' = exp[^CA:A"^], ci = ai (39) 

fc=l k=l 

and A ^ oo. The hydrodynamic chain satisfies the same generating function (j^ . 
where (sf. (US))) 

oo oo 

/X = J]6fcA^ = exp[- J]c_feA% co = -ln6o (40) 

k=0 k=0 

and A — s> 0. 

All the positive commuting flows in field variables Ck are 

n 

dtn'^k — ^ ^ QmCfc+ra-Tn,X) /c = 0, i 1, i 2, ... 

m=0 

where the generating function of conservation law densities for arbitrary positive com- 
muting flow is all negative flows are 

n 

dt_„-iCk = ^^^mCfc+m-n,X) /c = 0, i 1, i 2, ... , 

m=0 

where the generating function of conservation law densities for arbitrary negative com- 
muting flow is (j23j) . All these negative flows can be obtained from positive commuting 
flows (see above) by the reciprocal transformation (see Remark of previous section). For 
instance, the first negative flow is 

dt_,Ck = e'^d^Ck-i, k = 0, ±1, ±2, ... (41) 
Remark 4 Obviously, values Uk and can be expressed from analogues of / T^) and 

(see (Hg); 

OO OO 



— = 1 + ^^fifcA"'^ = exp[— y~^CfcA~^], ai = —Ci, A — oo. 

^ fc=l k=l 

oo oo 

- = Y}kX' = exp[Y^c_kX''], bo = e'^ A ^ 0. 

^ fc=0 A;=0 

Thus, hydrodynamic chain (jHTj) can be expressed via different moments (a^, bk), see, 
for instance Remark in section 6. 

Relationship (j39j) positive moments and positive conservation law densities can 
be expressed explicitly by next four recursive formulas, where the first of them 

k 

dok+i = y^^amdck+i-m + dck+1, k = 0,1, 2, ... 

m=l 
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is consequence of three local symmetry operators acting on space of conservation law 
densities a^: the shift operator 



5afc+i = = afc, /c = 0, 1,2, 

OCi 



the scaling operator 



mcm-^ — = kuk, fc = 0,1,2, 



m=l 



and projective operator 



°° d 
Sak = [ci + ^(m + l)cm+i^ — ]ak = {k + l)ak+i, k = 0,l, 2, ... 

m=l ^^"^ 

4 Finite-component reductions 

Theorem 5 T/ie hydrodynamic type system (0) w^/i arbitrary number of components 
is embedded into the hydrodynamic chain ^37\ ) if and only if 

N 

V. = f\{r')-c„ ci = 5^^„(0, (42) 

m=l 

where fi{r^) and ip j^{r'') are arbitrary functions. 

Proof. Any reductions are compatible with given hydrodynamic chain (jHTj) if every 
Cfc can be expressed as a function of just N independent Riemann invariants r*. Then one 
obtains 

VidiCk = diCk+i - CidiCk, i = 1,2, A^, = 0, ± 1, ± 2, ... 

It is easy to see, that 

diCk+i = {Vi + c.fd^ci, 1 = 1,2, N, k = 0, ±1, ±2, ... 

Thus, the second derivatives 

d,[{V, + c,)''^,c^] = d,[{V, + c,)''^,c^], i = I, 2, N, k = 0, ±1, ±2, ... 

yield the general reduction (ji^ . ■ 
Thus, all moments are 



N 

Ck 



E / [/™(A)]'-'rf^^™(A), = 0, ± 1, ± 2, ... 
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Remark 6 The hydrodynamic type systems 

rl_,=W\v)rl, 1 = 1,2, ...,N (43) 

embedded into first negative flow (see can he found in the same way (see ^^): 

from W'diCk = e^°diCk-i one can obtain diCk = (^)*'9jCo. In comparison with dtCk = 
{Vi + Ci)^diCo one can obtain 

N 

Next task is how to write hydrodynamic reductions in closed form via special variables 
like conservation law densities (see, e.g. ()30|) ). It means, that all higher moments a^+k 
must be expressed via lower moments (fc = 1, 2, ... , N). In particular case (fT^ . all 
higher moments a^+k are polynoms, which can be found from relations (see (|THl ) 

= [(1 + Aai + \^a2 + ... + A^a^v + A^+^a^+i + = 1, 2, ... 

For example, the first higher moment Qn+i can be found from more compact recursive 
relation 

r(i-i/.).r^) 
^r(i-fc-i/5)r(A:+i) ' 



(N+l) 

N+l-k-ik _ k I \u^k-l^ I I \N+l-k(N+l) 



where z^. is coefficient of series 



[a, + a^X + asA^ + ... + a;v+2-fcA''+'-1' = al + Xka^-'a^ + ... + z'^^^^> + ... 

All higher moments a^j^k are polynoms of lower moments a^, e.g.: 

a, = l±^a,a, - ^i±f^a?, N = 2, 
a, = '-^{2a,a, + al) - ^^^^^afa, + (i+^)(W(i+3.) ^4^ ^ ^ 3^ 

as = ^(aia4 + 0203) - ^^"^^^2"^^^'' (^1^3 + aia|) + 

(l+£)(l + 2£)(l+3£) 3-, _ (l+£)(l + 2£)(l+3£)(l+4£) 5 - 4 

and so on... First exceptional case, is chromatography phenomena {e = —1), then 

as = -aia2 + = 2, 

a4 = — i(2aia3 + a^) + \a\a2 — ^^af, = 3, 

«5 = -(0104 + 0203) + ^(0^03 + aial) - |afa2 + ^af, A^ = 4; 

the second exceptional case is (dispersionless limits of Coupled KdV and Coupled Harry 
Djth) e = —1/2, then first higher moment ak+i is just quadratic expression via lower 
moments a^. It means, that corresponding hydrodynamic type systems have at least one 
local Hamiltonian structure (see actually, the most number of local Hamiltonian 

structures is (A^ + 1), iff e = —1/2, see [IH]))- The third exceptional case is £ = — 1/M, 
where M = 3, 4, ... Then all higher moments will be quickly trancated (see above). 
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5 Commuting flows and reductions 

The corresponding linear system for conservation law densities (see jSEI) is 

^^'^^ = f( t'K\ k^ ^^^ - f( t'^^'lf kA h, I ^ k. (45) 

The general solution of such system is determined by N functions of a single variable (also, 
see |Sni)- In a particular, but very important, case of "e— systems" (see, for example, (fT^ ) 
this linear system is exactly A^— component generalization of the Euler-Darboux- 
Poisson system (see |22]) 

dikh = —^-^[dih - dkh], i 7^ k. (46) 

At first, it is necessary to find a generating function of conservation law densities /i, which 
can be found in comparison with (j^^ and ()42|) 

Mr,A)=expf (47) 



k: 



This formula (j47|) simplifies in case of "e— systems" (see (|14j) and (|15p. that was well 
known in case N = 2, for example, see [7j). Velocities of the commuting flows (i.e. 
Riemann invariants r* are considered as the functions simultaneously of three independent 
variables x, t, r) 

ri = w\v)rl z = l, 2, AT (48) 
can be found as the solutions of another linear system (also see |36j ) 

'""'^-A(?rW)<™' '''' 

Theorem 7 Any solutions of the linear system ^y{ ) are connected with the solutions of 
another linear system (sf. \4^ ) 

-"/.('•■) -A(r')^-''"'/.(r-)-A(r')*''- ' ^ " 

by the differential substitution of the first order 

w' = \fdih. 

Thus, the generating function of solutions (for the hydrodynamic type systems (jHI) and 
(1121), and (jH}; cf. (|T!^ and (jSSl)) by Tsarev generalized hodograph method (also see 
is 

TV AT r™ 

X + [/.(rO - Y.Mr^)\tr + TIF) ^^P[ S / = 

fc=l m=l 

(50) 

1 



irexp 



k=l^ A-/.(A) 
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where the generating function of commuting flows (see (|19p) in Riemann invariants is 

1 



^^'^ (l-/.(rO/A)/i 
Remark 8 // one selects monoms 

where Sk are arbitrary constants, then the generalized "e— systems' 

N 

r\ = y - Ve^r'^lrl. 



J2^mnrl (51) 

m=l 

has the generating function of conservation law densities (cf. and 

N 

^^=X[{l-rJ\)-^-. 



m=l 



In general case (N is arbitrary) the hydrodynamic type system i51]} is the natural N— parametric 
reduction of the hydrodynamic chain (|J7| ). When N = 2 this system l\51]} is natural two- 
parametric generalization of gas dynamics (0)- If we choose (cf. ^EE)) 



N , N 



Co = X^£™lnr^ cfc = i^£^(r™)^ k = ±1, ±2,... 



m=l m=l 



then the hydrodynamic type system satisfies for hydrodynamic chain |5 

Remark 9 In the particular case fiir^) = Si (Si are arbitrary constants and Si ^ Ek for 
i 7^ k; this is N— component generalization of gas dynamics, when the adiabatic index 
J = I) hydrodynamic type system 

N 
m=l 

is "trivial" (also ipi{r^) 7^ const^. In this case, the linear system has constant coeffi- 
cients 

d^h _ 1 dh dh 

where = ipi{r^). 
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6 General solution 



Description of the general solution for the linear system ()45p is a very complicated task. 
Construction of the general solution has been made only in case of "e— systems" (see j29j). 
when is arbitrary; cases N = 2 and N = 3 were completely investigated by G. Darboux, 
L.P. Eisenhart and T.H. Gronwall (see [7j and ^Hl)- The basic idea of how to construct 
a general solution (parameterized by N functions of a single variable, see ^|) of any 
over-determined linear systems like was presented in jSH] by recursive application 
of symmetry operators compatible with such systems. However, here we establish an 
alternative approach in the spirit of G. Darboux (see IHEI; also, see the section concerning 
elliptic coordinates in [7j). Elliptic coordinates /i„ (a = 1, 2, A^) appear in the theory 
of integrable hydrodynamic type systems associated with hyperelliptic curves, i.e. with 
"e— systems", where e = —1/2. G. Darboux suggested to introduce special variables r'^ 
{k = 1, 2, A^) for separation of coordinates in Laplace equation by the following rule 
(see (USD, when e = -1/2) 



N 

n (7a - 
2 fc=l 
/^a = 



n (7a - 7/3) 



N 

n (r' - 7/3) 

/3=1 



where 7^ are arbitrary constants (the denominator Yl {la ~ Ip) the first expression 

is just a constant multiplier, which does not affect the property to be a conservation law 
density, also see (j^^ ) and flat (not constant) metric 5'**(r) determined by 



TV 



N 



k\2 



0=1 



k=l 



Thus, elliptic coordinates coincide with the Riemann invariants for "e— systems", where 
e = —1/2. It is easy to generalize Darboux coordinates /i„ to arbitrary e 



N 



n (7a - r'' 

(/^a)-^/^ - ^ 



n (7a - Ip) ' 

P^a 



a = 1, 2, 



N. 



In this case, all "e— systems", for instance, (fT^ and (jHHj) can be written explicitly via /i^ 
in the conservative form (see (j30|l and p5|l in the particular case e = —1/2) 



'' N ^ \ 

^^M'^^' + 7a - £^lp ) f^c 
\ P=l P=l J 



(52) 



N 

Uilp) 
p=l 

7a 



-dr. 




^a 



a = 1,2, 



N. 
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Remark 10 Hydrodynamic type systems M'J\) and \3'J\) for another set of moments = 

N 

"^Zil pYil-'-p)"^^'^ (^^^ ^SM) cffT- be written as the following hydrodynamic chains 

13=1 

N 

dt.Ek = d,Ek+i + e[Eo - J^JplEk,, - E^Eo,^, k = 0, ±1, ±2, ... , (53) 

/3=1 

N 

dt_,Ek = l[{'yl,r[{l-E_^y^,Ek-l + {l-E_^r-'Ek-lE_,J, (54) 

/3=1 

where 

N N 

/3=1 /3=1 

Remark 11 The hydrodynamic chain i5^) is the same as the hydrodynamic chain HyU^) . 
because these two chains are connected by invertible transformation (see ^EB) cind ^EW) 

oo 

f,-'/^ = l + Y,EkX-'^'^'\ (55) 

where A — oo. The hydrodynamic chain ^54\) is the same as the hydrodynamic chain i35]} . 
because these two chains are connected by invertible transformation (see and 

oo 
k=l 

where A — > 0. 

Thus, our approach is the following: we mark arbitrary points A = 7q, (A^ distinct 
punctures) on the Riemann surface -F(A, /i) = (see (jlTj) and cf. (O); then we obtain 
special set of coordinates 

/^. = expf/^^, (56) 

which in fact is fundamental basic of linearly independent solutions for the corresponding 
linear system It means that any solution of linear system ()45|) can be presented as 
a linear combination of the basic solutions ()56|) with some coefficients. Finally, we just 
mention, that any over-determined linear system like (^3)) must have a general solution 
which depends on N arbitrary parameters 7^,. In our case, we should take infinite 
serieses of the conservation law densities /i^,^ {k = 1, 2, ...) starting near already fixed 
punctures 7^, 



=/ic. + (A-7„)/i„i +(A-7„)V2 + --, A^7„, a = 1, 2, ... , iV. 
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Thus, the general solution of linear system is 



N ' 



where fjsiX) are arbitrary functions, and the general solution of the hydrodynamic type 
system is given in an implicit form (see (jHl), (jl2I), P^jl - and (jnOj)): 



N , ^ 



N 



m=l 



/m(A) 



1 B 



E/^,(A)/i(^)(A)dA 



(57) 



where 



exp 



^^fc(A) 
7a - /fc(A) 



/x(r. A) = exp 



N 



c^^fc(A) 
A - /,(A) 



/ia + (A - 7a)/ia,l + (A - 7a) /^a,2 + 



A^7a, a = 1,2, 



N. 



The general solution of a linear system (jl^ can be presented in most possible explicit 
form in special case, when values e are integers for "e— systems" . The case N = 2 (namely 
Euler-Darboux-Poisson equation) was completely investigated (see, for instance, ISi])- Its 
generalization on A^— component case (|46|) . or moreover on case of arbitrary integers e.m 
(see dnH)) 

dikh = ziekdih - Sidkh], i ^ k. (58) 

can be made in the same way as in j^S]. For simplicity, here we shall restrict our consid- 
eration on case of (jlUjl (see [2U]). 

The general solution of (if e = ±n, n = 1, 2, ...) is 



N 



l{n) 



k=l ^ 



k\n 



m^k 



N 



k=l 



N 



m=l 



where (pi^{r^) are N arbitrary functions of a single variable (if we replace v^fc(A) 
V'i"^^^''(A) in second negative case, then all integrals can be expressed via finite num- 
ber of derivatives only). Thus, indeed, these solutions are general (see (2^) for positive 
and negative integers e (see flT^ and (jlHl))- The general solutions for can be obtained 
by recursive application of Laplace transformations (see |llj) to above formulas (also, see 

IB) 
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If e is negative and e 7^ — n, n = 1, 2, ... , then above-mentioned solution (right) 
easily to generalize 



N ■'"'1 N 



k=l-^ m=l 

"Ik 



where 7^, {k = 1, 2, ... , A^) are arbitrary constants. If e is positive and e n, n = 1, 
2, ... , then above-mentioned solution (left) easily to generalize just in case, when eN is 
integer, then 

'-'q n (A- r»')' 

m=l 

where (fc = 1, 2, ... , A^) are simple small contours surrounding the points X = r'^ 
{k = 1, 2, ... , A^). However, in general case (when e is positive and eN is not integer) 
these contours could not be closed on corresponding Riemann surface, because a sum 
of all phase shifts (for every point) will not be proportional to 2itM, where M is some 
integer. For avoiding this problem one can introduce another set of contours-dumb-bell 
shaped figures C^, k+i surrounding every two neighbor points X = r'^ and A = r^"*"^ (A; = 1, 
2, ... , N). So, integration must change sign twice from clockwise to anticlockwise, then 
every time phase shift will be 2n exactly. However, the number of contours must be equal 
to A^— 1, because in opposite case (if number is A^) all contours became hnearly dependent. 
Thus, in this general case a general solution of ()46|1 parametrized by A^ arbitrary functions 
of a single variable is 



7V-1 







N N 
m=l m=l 



where for simplicity we assume (without lost of generacy), that real parts of Riemann 
invariants (branch points on a complex Riemann surface) are numerated as follows: 



< Rer^ < Rer^ < ... < Rer^. 



7 (2+l)-integrable hydrodynamic type systems 

Benney momentum chain is equivalent to the hierarchy of (2 + 1) hydrodynamic type 
systems embedded in dispersionless KP hierarchy as Khohlov-Zabolotskaya equation 

{ut2 - uux)x = ut^ti, (59) 

which can be obtained from the coupled equations of the Benney momentum chain 
and one equation of its first nontrivial commuting flow (see the next section) 

dt,Ak = dxAk+2 + AoAk,x + {k + l)AkAo,x + kAk-iAi^x, /c = 0, 1, 2, ... (60) 
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by eliminating moments Ai and A2: 

dt.Ao = d^A,, dt,A, = d,[A2 + ^Al], dt.Ao = ^,[^2 + Al], 
where u = Aq. 

Let us start now with the hydrodynamic chains (jHUj) and (jH^ . ehminate field variable 
03 from couple equations from first hydrodynamic chain (|3U|) and one equation from the 
second hydrodynamic chain (j6U|l 

dt^ai = d^[a2 - al], dt^a2 = d^la^ - aia2\, dt^ai = d^la^ - 2aia2 + al]. 

Then we came to a new integrable (2 + 1) hydrodynamic type system 

Uti = Wx, Ut2 = Wti + UWx - WUx, (61) 

where 

u = ai, w = a2 — a^. 

It is easy to check, that all possible hydrodynamic reductions of this system ()6H) (see, for 
instance, approach in [25] ) 

rl=f^,ir)rl, rj^ = C,(r)r;, z = 1, 2, ... N 

are completely the same as those found already (see (jl^ ). where 

Q. = f-{rl - /.(r^)E^™(r™) + - ^^^^(r") " E / ^(^)#.(^)- 

fc=0 \fc=0 / k=0'^ 

Moreover, one can obtain a whole hierarchy of such integrable (2 + 1) hydrodynamic type 
systems like ()61|) by eliminating some other field variables in combination with another 
commuting fiows of hydrodynamic chain (jHU)) . For example, two another equations (see 
dnZI) and dUD) 

5tie-^° = 9,[-Cie-^»], dt_,c, = 9,e^° (62) 

yield a new integrable (2 + 1) hydrodynamic type system (its (1 + 1) hydrodynamic 
reductions are exactly). 

8 Another hydrodynamic chain 

Now we start with the integrable hydrodynamic type system [27j 

^ = 1, 2, ... , iV, (63) 

when Ek are arbitrary constants. This system can be rewritten as the hydrodynamic chain 
dtCk = cid:,:Ck - codxCk+i, k = 0, ±1, ±2,..., (64) 

where moments 

N 

Ck = 5^r'"(£„)^ 

m=l 



N 

E 

m=l 



Si 



N 
r?i=l 
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Theorem 12 Under the reciprocal transformation 

dz = — dx + —dt, dy = dt 

Co Co 

this hydrodynamic chain linearizes 

dyCk + d^Ck+i = 0, k = 0, ±1, ±2,... (65) 

It means that any reductions such as of the hydrodynamic chain (jMj) linearizes 
under above reciprocal transformation. The solution of the hydrodynamic chain is a 
set of the separate Riemann-Monge-Hopf equations 

rl + Mryi = 0, 1 = 1,2,... ,N, 

where /j(r*) are arbitrary functions. Thus, every integrable reduction of hydrodynamic 
chain (jMjl has the simple form 

ri = [cofi{r') - ci]ri, 

where 

N 

Cfc = E / [/™(A)]'t/^„(A), k = 0, ±1, ±2, ... 

and iprni^"^) arbitrary functions (by scaling tlj^{r"^) Rm integrable hydrodynamic 
reductions are parameterized by arbitrary functions of a single variable only). 



9 Open Problems 

The Benney momentum chain (the Zakharov reduction) is a dispersionless limit of the 
vector nonlinear Shrodinger equation (see (jlll) . (fT^ and [HHI)- The inverse problem is: 
how to reconstruct a dispersive integrable analogue of given hydrodynamic type system. 
A dispersive analogue is known (Coupled KdV is a dispersive analogue of system (fT^: 
Couple Harry Dym is a dispersive analogue of system (IHHj). see, i.g. |^ and [TB]) just 
in case of "e— systems" with e = —1/2. The KP hierarchy is a dispersive analogue for 
the whole Benney momentum chain (Q) (as KP equation is a dispersive analogue of the 
Khohlov-Zabolotzkaya system (j59j) ). but similar dispersive (2 + 1) analogues for the whole 
hydrodynamic chain ()37|) or for (2 + 1) hydrodynamic type systems (1611) or ()62p still are 
unknown. 

Local Hamiltonian structures for the hydrodjTiamic type system p3|) were completely 
investigated in [221 and ^H]- It was proved, that if = 2, then the hydrodynamic type 
system (jl3p for any e has three local Hamiltonian structures (also, see j28j and ^3)! 
A^ = 3 and e = —1/2, then it has four local Hamiltonian structures; if A^ = 3 and s = 1, 
then it has two-parametric family of local Hamiltonian structures (also see jHS] and [T^l: 
if A^ > 3, then e = —1/2 and it has (A^ + 1)- local Hamiltonian structures. 
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Hamiltonian structures of integrable hydrodynamic type systems are determined by a 
metric gu (see, [9J. The metric 



/ = Ci('^*)exp 



■2E 



with an arbitrary functions Cj('"*) determine a nonlocal Hamiltonian formalism (see 
and [36]) of hydrodynamic type systems (jH)), P2|l . P3|) . (jSj) and their commuting flows. 
Unfortunately, local and nonlocal Hamiltonian formalism has been done just for the hy- 
drodynamic type systems ()13|1 when e = ±1 and e = —1/2 (see for instance DEj and 
|12j). However, the problem of a description of local and nonlocal Hamiltonian structures 
in general case ()42|) still is open. Nevertheless, this problem can be solved by the Dirac 
restriction of a Hamiltonian structure (see for the beginning jT^) known for the whole 
hydrodynamic chain, as it was already done in jB] for another hydrodynamic chains. 

Starting point of such investigation is a Lax-like representation. For instance, the 
Lax-like representation (see Q) for the dispersionless KP hierarchy (i.e. the Benney 
momentum chain (^) is well known (see 

d,.X = {Qn, A} = ^- - n = 0, 1, 2, ... , (66) 

where Qn is the part, polynomial in /i, of A". Also, the first local Hamiltonian structure 
for whole Benney momentum chain 

dt^Ak = [fe^fc+m~l^x ^f?~^^ + ("^^fc+m-l^^^)x], (67) 

where the Hamiltonian is if2 = | /[^2 + AQ]dx, was constructed in [23| (the relationship 
between formulas and (jU7|) was found in too; first notrivial commuting flow 
(see fl60p ) is determined by the next Hamiltonian = ^ J [^3 + 3y4ov4i](ix; functional 
Hq = J Aodx is a Casimir of this Hamiltonian structure, the functional Hi = J Aidx is a 
momentum of this Hamiltonian structure). 

Similar Lax-like representation for the hydrodynamic chain ()37p was established in ^ 

(cf. dnnD) 

dL do 

dt^L = (Q„, L) = g„— - ^L, n = 0, 1, 2, ... , (68) 

where 

Qn = (A"L)+, L = l + Go/A + Gi/A^ + Ga/A^ + ... 

The corresponding first hydrodynamic chain 

dt^Gk = dxGk+i + GoGk,x — GkGo,x, k = 0, 1, 2, ... 

/ N \ 

is exactly the hydrodynamic chain (jK!^ . where e = 1, and linear term — e Yl'la ] -^k x 

V/3=i J 

f N \ 

is removed by shift of independent variable {x ^ x — e \ J2li3 ] '^)- Thus, a generating 
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function of these moments is (see (f^. and (|HHj)) 



oo 



k=0 

The ahernative Lax-hke representations are 

k k 

dt,P = ES^A^-'"^,, p], dt_,_,p = [J2bmX"'~'''d,, p], k = 0,l, 2, 



m=0 771=0 



where p = p^ (see 

We suppose that the hydrodynamic chain (jHTj) and its commuting flows have local 
Hamiltonian structure (sf. (jU7j) ) 

k=l 

where are some functions. 

The Hamiltonian structures of integrable hydrodynamic type systems can be success- 
fully investigated by apphcation of methods from the differential (see, for example, 
j2ni-[3Hl, [El, [ISl, [12] and [in]) algebraic geometry (see, for instance, [Hj and [22|)- An 
alternative way is following: assume that our given integrable hydrodynamic type system 
{N components) is a some reduction of some "bigger" integrable hydrodynamic type sys- 
tem {N + M components); assume that Hamiltonian structure of such "bigger" integrable 
hydrodynamic type system is already known. Then the direct application of the Dirac 
restriction to this Hamiltonian structure (the first step in such procedure (see [12]) is the 
choice of some Riemann invariants r'^ = const, k = 1, 2, ... , M) yields the transformed 
Hamiltonian structure of a "restricted" hydrodynamic type system. The Dirac restric- 
tion of Hamiltonian structures (in algebraic language) was developed in application to 
hydrodynamic chains and their reductions (see [Uj). The first step in such procedure is 
the recalculation of the Lax-like representation (such or (j68j) ) to the Hamiltonian 
structure of a whole hydrodynamic chain (see for instance ()67|)). 



10 Conclusion 

In this article we present a recipe: how to construct a hydrodynamic chain starting from 
any given hydrodynamic type system with polynomial (or rational) velocities with respect 
to their field variables (for simplicity we have mentioned just two cases: namely Benney 
momentum chain, whose moments are connected directly with some conservation law 
densities {uk, rj^.) and the hydrodynamic chain (jHTj), whose moments are connected directly 
with Riemann invariants) . In fact, it means that any integrable hydrodynamic type 
system (written in Riemann invariants) with such polynomial velocities is 
embedded in hydrodynamic chain (|n7j) or its higher (or lower) commuting 
flows. Any integrable hydrodynamic type system written in Riemann invariants with 
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rational velocities 



(r.)"+g.(r-)"-' + ...+gM ^,. ^ _ J 2 ^ 

"^i 9»(r<)A- + ei(r')f-i + ... +eK " LA-,™ (DJ) 

has a generating function of conservation laws 

^ r A +giX + ...+gM . 
A + CiA + ... + ex 

For simplicity we assume that coefficients and of such rational velocities are some 
symmetric (not necessary to be polynomialsl) functions of Riemann invariants; K, M 
and are arbitrary natural numbers. As example, we can take any integrable systems 
embedded into 2x2 spectral transform like Korteweg de Vries equation, Bonnet equation 
(Sin-Gordon equation) and nonlinear Shrodinger equation. All their Whitham deforma- 
tions (i.e. hydrodynamic type systems, see, for instance, [HI, |22]) have such representation 
(jnni) as consequence (fTHjl written in abelian differentials on hyperelliptic surfaces 

dtdp = dj:dq, 

where 

dp = 1 dX, dq = go dA, 

N N 

n (r^ - A) JU (rm - A) 

m=l V m=l 



and velocities of (jU^ are 



jr')^ + gi{r')^-^ + ... + gM ^ dq 
^0 (^ri)K + ei(r^)^-i + ... + ~ dp '^="' 



Substituting (sf. (jSl) 



yU = 1 + y~]afcA 



k=l 

into generating function (fTOj) one can obtain similar formulas and results as it was done in 
this paper. Next step is replication of integrable hydrodynamic type systems as different 
hydrodynamic reductions of hydrodynamic chains. The main advantage of such repli- 
cated systems is preservation of some properties of original hydrodynamic systems like 
generating functions of conservation laws and commuting flows (see (0)) and (jHI) for 
Benney momentum chain see (fT^ . (j^^ and (j^Uj) for the hydrodynamic chain (jHTj) ). 
Thus, a problem of integrability is much simpler - all that necessary to do: to construct 
a general solution (starting from already obtained generating function, see ((Tj), (pH) and 
(jnZj)), parameterized by A^ functions of a single variable (see [Sni) and to solve Cauchy 
problem, that in fact is done to this moment just for four cases (the Zakharov reduction 
of the Benney momentum chain, see ^7j; linearly degenerate system, see j^ZI; this par- 
ticular class is "e— systems" with e = 1; hydrodynamic type systems of the Tample class. 



23 



this particular class is "e— systems" with e = —1; Whitham hydrodynamic type systems 
related with hyperelliptic surfaces such averaged A^— phase solutions of Korteweg de Vries 
equation (KdV) or nonlinear Shrodinger equation (NLS), it was done in articles of G. El, 
T. Grava, B.A. Dubrovin, F.R. Tian, J. Gibbons and many others). 
Moreover, any two commuting flows of (jSHI), for example (PT|) 

k n 

dz = fi[dx + ^dm\''~'"'dt^ + ^S„A"-'"dt"], k,n = l, 2, ... 

m=0 m=0 

yield hydrodynamic type systems with rational velocities 

k 



tk ~ 



i\k—m 



rl = ^ Mn, z = 1, 2, AT. 



m=0 

A more complicated rational dependence can be obtained by application of a generalized 
reciprocal transformation (see, for instance, [12] and ^H]), starting from (j^^ and its 
commuting flows. 

Thus, this is powerful tool for classification of integrable hydrodynamic type systems 
and their integrability. Moreover, if any given hydrodynamic type system has a gen- 
erating function of conservation laws (see for instance, ^ or (j^^ ). it means that the 
corresponding hydrodynamic chain has the same generating function. For example, if a 
some hydrodynamic type system has the same generating function as the Benney mo- 
mentum chain (0]), it means that this hydrodynamic type system is a some reduction 
of the Benney momentum chain. Thus, this is a wonderful symptom in recognition of 
an immersion of any unknown hydrodynamic type systems into already known hydrody- 
namic chains. Thus, if one can construct a generating function of conservation laws for 
some hydrodynamic type system, it means that, in fact, hydrodynamic chain is already 
constructed (and may be recognized, because, obviously, amount of hydrodynamic chains 
is much smaller than amount of integrable hydrodynamic type systems). 

However, the problem of a description of all possible reductions is very complicated. 
For instance, this problem for the Benney momentum chain is still open (see PH])- How- 
ever, this problem for the Boyer-Finley momentum chain (continuum limit of the Darboux- 
Laplace chain, which also is known as two-dimensional Toda lattice, see [23 and |22|) 
in fact is not exist, because both mentioned hydrodynamic chains are related by special 
exchange of independent variables, see |21j). Thus, we are lucky to solve this problem for 
the hydrodynamic chain ^7\i . 

Acknowledgement 13 / am grateful to my friends and colleagues Eugeni Ferapontov for 
his suggestions in improvement of my calculations and Yuji Kodama for useful discussions. 
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